global C1 C2 C3 dif ky kz
Te = 150; % in eV
mi = 2*1.67e-24; % in g
me = 9.1e-28; % in g
q = 1.6e-19*3e9; % in statC
B = 20000; % in GS
c = 3e10; % in cm/s
n = 1e20*10^-6; % in cm^-3
lnA = 22.36+1.5*log(Te)-0.5*log(n); % Coulomb logarithm
h = 3; % boundary layer thickness in cm
kappa = 1/170*h; % normalized curvature 
chie = 3.2*1.6*3.44*10^-7*Te^2.5/me/n/lnA; % electron heat conductivity in cm^2/s
chie = chie*1e-5;
wci = 9.577e3*B/2; % ion gyrofrequency
ti = 2.085e7*sqrt(2)/lnA/n*(Te^1.5); % ion collision time
nu_i = 3/10*1.6e-12*Te/wci^2/ti/mi; % ion kinematic viscosity in cm^2/s
% nu_i = nu_i*10^2;
dif = 10; % Prandtl #
% dif = 10*dif;
% Ra = kappa*Te*1.6e-12*h^2/(mi*chie*nu_i);
% C2 = 0.1*C2;
t_norm = h^2/chie;
gradT = -1; % normalized temperature gradient (constant)
Ra = 20000;
C1 = Ra*dif;
C2 = dif;
ky = 3; % used in T0
kz = 0.005;